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Abstract 
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1 Introduction 



Quantum Chromodynamics (QCD) is formulated in terms of a functional integral for quarks and 
gluons. The observed particles, however, are mesons and baryons. For a computation of mesonic 
properties one may imagine a two-step procedure: First one integrates out the gluons and obtains 
an effective functional integral for the quarks. Second one computes the mesonic properties from 
the effective quark model. At first sight this seems an extremely difficult task since the result of 
the first step must be a complicated nonlocal effective action for the quarks. For example, the first 
step already contains all the information needed for a determination of the heavy quark potential 
in the quenched approximation. Nevertheless, the following four observations indicate that this 
task may not be insoluble: 

(i) Due to asymptotic freedom the momentum dependent multiquark interactions can be com- 
puted approximately for momenta sufficiently above typical strong interaction scales, e.g. the 
"confinement scale" Aqcd- The full understanding of non-perturbative gluon physics is only 
needed for the low-momentum behaviour of the quark n-point functions which are obtained in 
the first step. Furthermore, non-perturbative flow equations based on the scale evolution of the 
effective average action can extrapolate into a momentum region beyond the perturbative range. 
The effective average action F^. includes in a Euclidean formulation all quantum fluctuations with 
momenta p"^ larger than some infrared cutoff k"^. It is the generating functional for the IPl n-point 
functions. A solution for k ^ accounts for all quantum effects and therefore contains the rele- 
vant information for particle masses and interactions once momenta are analytically continued to 
Minkowski space. The dependence of the effective average action on the scale k is described by 
an exact renormalization group equation 0] . This allows us to explore the complicated infrared 
region gradually. 

(a) The functional integration for the quarks in the second step can also be done in small 
consecutive pieces by the use of exact renormalization group equations. It has been shown that a 
formulation in terms of nonperturbative flow equations derived from the effective average action 
0, 1^ can deal with nonlocal quark interactions 0. 
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Combining this with (i) one expects that a rehable computation of the effective average action 
for quarks (after integrating out the gluons completely) should be possible at least for momenta 
and k larger than the vicinity of the confinement scale. 

(in) A first investigation of the flow of the effective average action for quarks has revealed 0] 
that mesonic bound states typically appear at a "compositeness scale" k^p ~ 600 — 700 MeV. This 
scale is considerably above the confinement scale and the strong gauge coupling is not yet very 
large. For example, the loop expansion of the heavy quark potential still seems to be qualitatively 
valid for y/q^ ~ 700 MeV [^]. The prospects that at least the physically relevant parts of F^^ can 
be computed rather reliably seem therefore quite promising. This quantity is interesting in its 
own right since it may be considered as the effective action for Nambu-Jona-Lasinio type models 

i- 

(iv) It is possible to rewrite the functional integral such that at the scale k^p the composite 
mesonic degrees of freedom are included as explicit fields The phenomenon of spontaneous 
chiral symmetry breakdown occurs as k is lowered beyond /c^ (typically for k ~ 450 MeV) 0. 
As a consequence, the quarks acquire a constituent mass rriq ~ 300 — 350 MeV and decouple from 
the flow of the mesonic interactions as k decreases below ruq. The uncertainty in the strong non- 
perturbative confinement effects expected for the quark n-point functions with momenta g^^Ag^^^.^ 
affects only moderately the meson physics. Furthermore, the running of many mesonic couplings 
is governed by a fast attraction towards partial infrared fixed points [0]. Except for a few rele- 
vant parameters many properties of the meson physics become therefore rather insensitive to the 
precise form of F^^. 

This paper aims for a computation of the effective four quark interaction in the first step of 
such a program. We will be satisfied with quantitatively reliable estimates for momenta down 
to i/g^ ~ 300 — 500 MeV whereas larger uncertainties are tolerable for lower momenta. More 
precisely, we want to compute the effective average action for the quarks at a scale k^ (typically 
near the tau lepton or charm quark mass), with the gluon degrees of freedom completely integrated 
out and the quark fluctuations with momenta larger than k^ included. 
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Our main tool is the effective average action [Hi- It is a functional of the fields where the 
derivatives with respect to x ^-re the full or proper vertices. The average action Tk[x\ obtains by 
adding to the classical action an infrared cutoff piece which is quadratic in the fields 

Here Rkip"^) constitutes a momentum dependent infrared cutoff. It acts like a mass term i?^ oc k"^ 
for the low momentum modes and is exponentially suppressed for the high momentum modes 
{Rk oc p2g-p^/fc^ fQj, p2 ^ ^2^_ r^^^ effective average action is then computed by the usual steps 
(introducing sources and performing a Legendre transform) while keeping the infrared cutoff 0. 
As the infrared cutoff scale k is lowered, quantum fluctuations with larger and larger length scales 
are effectively included in the functional integration. The variation of F^, with k is described by 
an exact functional differential equation (t = Infc/fco) @ 

In fact, the only fc-dependence of F^ arises through the quadratic infrared cutoff (P. It is therefore 
proportional to dR^/dk and to the connected two point function (t^^^ + -Rfc) • One can use 
different cutoffs ksik) and kp{k) for the bosonic (gluons and ghosts) and fermionic (quarks) parts. 
Correspondingly, (^f[.^^ + Rk) denotes the projection of the full field dependent propagator on 
the bosonic subspace and similar for the fermions. (For details see reference [Q.) The traces in 
eq. (0) contain a momentum integration and a summation over internal and Lorentz indices. Our 
exact equation has the form of a renormalization group improved one loop equation. The only 

(2) 

difference from a one loop equation is the appearance of the full inverse propagator F^ , i.e. the 
matrix of second functional derivatives of F^ with respect to the bosonic and fermionic fields Xi 
instead of the classical inverse propagator S**^^^ . This turns equation (^ into a functional differential 
equation and renders it exact to all orders in the loop expansion and including all nonperturbative 
effects. Furthermore, the momentum integration in equation (0) is both ultraviolet and infrared 
finite. Only a narrow range in the loop momenta ~ k'^ contributes effectively. In consequence, 
a computation of dtVk only needs knowledge of the physics in a momentum range near k. 
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The computation of the fermionic average action at the scale can be achieved by starting 
at some high scale ks = kp = ko within the range of validity of perturbation theory and lowering 
simultaneously the bosonic and fermionic infrared cutoff (fc^ = kp = k) down to k^, using equation 
@. One then keeps kp = k^ fixed and removes the infrared cutoff for the bosons {ks = k) by 
following the evolution towards ks 0. For this purpose one has t = Inks/ko and the second 
term in eq. @ does not contribute. Since the infrared cutoff k^ for the quarks is still present 
all loops containing internal quark lines are suppressed for k <^ k^ - similar to the situation for 
heavy quarks for Euclidean momenta. The resulting effective action Tq^^^ [x] is still a functional of 
both bosonic and fermionic fields. In order to obtain the effective action involving only quarks ip 
one has to solve the bosonic field equations from ro,A:^. The solution expresses the gauge fields as 
functionals of the fermions. (The solution for the ghosts corresponds trivially to vanishing ghost 
fields.) Inserting this solution into ro,A,.^ yields the fermionic effective average action rfc^[-?/'] which 
depends only on the quark fields. We want to compute in this paper the terms in [■?/'] involving 
two or four powers of ip. The second step of solving the flow equation for [■?/'] for A;^ — is 
outside the scope of this paper and we refer the reader to [If. 

We also relate the effective four quark interaction to the heavy quark potential. In section 
6 we present an approximative estimate of the Fourier transform of the heavy quark potential 
V{q'^) in the "quenched approximation" where no light quarks are included {Nf = 0). An optimal 
renormalization scheme guarantees that Vlq"^) almost agrees with the two loop potential in the MS 
scheme for 3A^ < ^ 10Aj;jg. (Here = 315 MeV is the appropriately adapted two loop 
confinement scale.) For lower momenta \fq^ < 2.5A(^ our result for the potential strongly deviates 
from the two loop potential. It turns out to be close to phenomenologically motivated potentials 
as the Richardson potential in the momentum range A^ ^ y/q^ ^ 2.5Aj^. Quarkonia level can 
therefore be described successfully with our approximate solution of the flow equations. We also 
attempt a first estimate of the string tension. The order of magnitude comes out satisfactory. The 
detailed quantitative value stronly depends on the truncation, indicating an insufficiency of the 
present approximations for very low values of g^. 
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Figure 1: The flow equations for tlie gluon (wiggly lines) and ghost (broken lines) inverse 
connected two-point functions r^.^-*. Here dt = ^^g^ denotes a derivative operator that only 
acts on the explicit fc-dependence in the cutoff-term. The full dots correspond to the full k- 
dependent propagators including the cutoff-term whereas the empty dots stand for the full IPl 
n-point functions. Combinatorical factors are not displayed. 

In the course of computing ro,A:^[x] we have to determine the propagators for the gluons and 
ghosts. Estimates for these propagators have been given using Schwinger-Dyson equations (SDE) 
1^. In a certain sense our non-perturbative flow equations can be interpreted as a differential 
form of the SDE. In contrast to the SDE only a narrow momentum window ~ k"^ contributes 
effectively to the momentum integration in the flow equation and truncations should therefore 
be more reliable. For a possible comparison with other approaches we also present in this paper an 
estimate of the gluon and ghost propagators in Landau gauge and the quenched approximation. 
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Figure 2: Momentum dependence of q^V{q^) for different values of tlie cutoff k. The dotted line 
corresponds to A; = 500 MeV, the dashed line is for k = 100 MeV, and the solid line displays the 
final result for k = 0. 



In practice, the most important ingredient of our calculation is the numerical solution of non- 
perturbative flow equations for the gluon and ghost propagators. These equations obtain from the 
exact equation (H) by a truncation which is explained in detail in the next section. Graphically, 
the exact system of flow equations for the gluon and ghost propagator (for fixed kp = k^) is 
represented in figure 1. The truncation corresponds to an approximation for the vertices and the 



precise form of the flow equation can be found in section 3 (equations (13) and (14)). As the 
infrared cutoff k = ks is lowered, more and more of the fluctuations are included. This results 
in a transition from an almost classical form of the potential V oc 1/g^ for k = 500 MeV to a 
non-perturbative form for — 0, as can be seen in figure 2. 

In order to connect our results for small with the perturbative results for large it is 
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important to choose the initial values for the numerical solution of the flow equation carefully. 
Otherwise, the overall (non-perturbative) mass scale appearing in the potential for the range 
300 MeV < i/g^ < 700 MeV would remain a free parameter. We discuss in detail the choice of 
initial conditions in sections 4 and 5. This allows us to establish an explicit relationship of the 
relevant mass scale with the perturbative mass scale appearing in the MS scheme. 

We should point out that the present paper is closely related to earlier work in the context of 



exact renormalization group equations on QCD propagators and the heavy quark potential |T0 
even though the present perspective is somewhat different. 

2 Effective four-quark interaction 

In this work we compute the effective four quark interaction in Fo.a;^. We always assume here 
kp = kjjj held fixed and identify k = ks- In order to approximately solve the exact fiow-equation 
@ we have to resort to a truncation of the most general form of the effective average action. We 
work in the general setting of flow equations for gauge theories [|1T], |T^, |13|, with Landau gauge 
fixing. This is motivated by the fact that a vanishing gauge fixing parameter a = is a fixed 
point of the flow equations |15| and will turn out to be crucial later. We work with vanishing 
"background fields" where the effective average action is not a gauge invariant functional. As 



far as the Slavnov- Taylor identities are concerned one notes that Tk contains counterterms, 
Tct, which are induced by the breaking of BRS-invariance due to the infrared cutoff Rk- For 
k^ikp — > the BRS invariance is restored and the counterterms vanish. We may implicitly define 
Vet by requiring that — obeys the standard Slavnov- Taylor identities. An exact solution of 
these identities which is sufficiently general for our purpose can be written in the form (a — » 0) 

Here Fj^^ is at this point an arbitrary gauge invariant functional of the rescaled gauge field A 
and the quark fields ip, with covariant derivative D^{A)ip = d^ip — igTzA^^ip. The fields A and 
^ are related to the gauge field A and the antighost ^ by a momentum dependent wave-function 
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renormalization 

A = H{-d^)A , l=H-\-d^)l (4) 
Comparison with (|) relates H to the inverse ghost propagator Ggh by 

G,n{q')=q'H-\q') (5) 

The function H[q'^) will depend on k and the scaling of A is such that the covariant derivative 
always involves a fixed gauge coupling g. This brings the present formulation close to the back- 
ground field formulation |]ll|, even though we will work here with a vanishing background field. 



From the ghost sector in our truncation we only need information about the two-point function 
Gghi.'f)) since the A^^-vertex is fixed by the constant g in Landau-gauge. 

In this setting the effective quark interactions can be extracted from the gauge invariant kernel 
rj,„j;[y4, ^/^j. For the general discussion we will consider the truncation^ 

a 



a,b y 

which will be further restricted below for our practical calculation. Here Zp, Yp, Z^, Y^, and E"-'^ 
are arbitrary functions of the covariant Laplacian (involving A) in the appropriate representation 
and F^j^ = dfj,Al — dyA^^ + gf^^^A^A'^. The index a labels the different quark-fiavors. We denote 
the contribution from T^ny to the transversal two-point function for A (not A) by 6*^(9^), 

GAiq') = {q'Zp{q') + q'Yp{q')) H\q') = {Zp{q') + q^Yp{q')) q'Gj{q') (7) 

The last term in equation oc E accounts for part of the IPl four quark Green-function. 

In order to extract the information on effective quark-interactions we have to solve the field- 
equation for the gluon-field as a functional of the quark field and insert the solution into the 



"'^We work in Euclidean space and 7 is the Euclidean analogue of 75 (for conventions see Ef 
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effective action. Restricting ourselves to the four-quark interaction, the solution of the field- 
equation for A derived from (|^) is needed only in lowest order in an expansion in {ip'il') 0- Upon 
inserting the solution into rj„i,[A, ■?/'] we find Teff[ip] a^ 



effi 



X 



X 



{g-d'Y4-d')Y 



-\ ab 



7b 



+ E{-d') 



X v.,7'^r^zJ(-9>, +... (8) 



where we have neglected terms with different index-structures for the four-quark interaction and 
higher powers of We will concentrate on the quantity 
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(Z^(g2)+g2y^(g2))^2 



where the normalization is such that V^(g^) coincides with the Fourier transform of the heavy 

quark potential^ V^(g^) in the one loop approximation up to terms oc g'^ /q^- (In equation (|^) we 

have suppressed the flavor indices of and E). We want to compute V{q^) in the limit A; — > 

such that all gluon-fluctuations are integrated out. We will also henceforth approximate^ the 

contributions from Y^ and E in the Landau gauge by a ^-independent constant A oc g'^. This is 

suggested by the observation that the running of Y.^ and E involves loops with internal fermion 

lines and is suppressed for a fixed fermionic cutoff kp- In summary, our calculation amounts to 

the determination of the scale dependent gluon and ghost propagators GA^q^) and Ggh{q^)- 

^ The additional ^JoV'a 4- vertex connected to a non-constant D^) contributes to the field-equation for A, 

but leads to a different index structure for the four-quark interaction. This also holds for gauge invariants of the 

type Tpal'^l'^ F^^TziJja which are not contained in our truncation. Nevertheless, there are in general also neglected 

invariants contributing to the index structure (V'aT'^V'a)^- 

^Here V{q'^) corresponds to the quenched approximation without light quarks {Nj — 0). 

^The solution of the modified Slavnov- Taylor Identities in ref. | [Tc| ] is not general enough to admit nonvanishing 
E or Y^. Note that in general gauges the expression for V does not reproduce the correct running of the gauge 
coupling even to one-loop. This is due to contributions that are absent in Landau gauge. 
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For > we also have to specify our truncation for the counterterms contained in F^. In 
principle, this piece is completely fixed by the additional terms appearing in the Slavnov- Taylor 
identities for k > [|I^. We only include explicitly^ a gluon mass term, given by 



Z,,Rf\x) 1 Z,,R^t\x)± {G,n{x) + ^ 



X 

n, ^ ax \ ^'-v/ V 

^ ZgHRi^\x) 

{Ga{x) + r\^\x)) (G,n{x) + R^f\x) 

1^ Z,nRi^\x)j-^{G,,{x) + R^t\<, 
~^\GA{x)+Rt\x))[G,,{x) + Rf^' 



[^x 

Here the wave function renormalization constants are given by 



y _ 9Ggh{q 

,2=0 '^y 



2^ 



f=0 



In the flow of general models for vector fields this mass term would be a relevant parameter, 
whereas in gauge theories its value is dictated by symmetry. Computing the mass term from the 
flow of truncated evolution equations could easily lead to large deviations from the "gauge theory 
hypersurface" in parameter space. It is therefore preferable to determine the mass term directly 
from the appropriate symmetry relations as in equation ([10|). As the infrared cutoff is removed 
the mass term vanishes for A; ^ 0. The mass term depends on the specific form of the infrared 
cutoff-function Rk- A natural choice for Rk{p^) is 

R^^p") = (12) 



^The contribution from Vet to the gluon two-point function Ga('?^) also includes q^-dependent terms. The flow- 
equation for Gyii((7^) presented below actually describes the evolution of the full gluon two-point function derived 
from Tim, + Tct, with the constant for zero momentum (the massterm m\) subtracted such that Gyi(O) = 0. The 
final result for Ga('?^) for k = Q nevertheless describes only Tinv, since the contribution from Vet vanishes. For the 
approximation of the vertices used in this work there is no need to separate the contributions from Tinv and Yd- 
Thus the g^-dependent part in Tct is implicitly included. Only the mass term needs to be determined in addition 
to GA{q^), since by definition GA{q^) — GA{q^) — Ga(0) ~ GA{q^) Similar remarks apply to Ggh where we 

note that Ggh{0) remains identically zero (cf. also eq. (^^) such that m'^h — for all k. 
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with different fc-dependent wave function renormalization functions Z for gluons and ghosts. In the 
region of approximately perturbative behaviour for k > 3Aj^ (with the two loop confinement 
scale in the MS scheme) there is also approximate scaling in the renormalized quantities and we 
choose fc-dependent Z in by identifying the constants in the infrared cutoff with the running 



wave function renormalizations (11) for the gauge fields and ghosts, respectively. For k < 3A^ 
scale invariance is strongly violated anyhow and we do not change Z with k anymore, i.e. we 
take Z = Z{3Ajj^) in Rk. For low momenta <^ k"^ the infared cutoff Rk acts like a mass 
term oc Zk"^. Typically, the size of the mass term contained in Td is smaller than Zk'^. In 
a certain sense, the mass term oc k^ actys similar to a modification of the infrared cutoff Rk- 
Including the contributions from Fct the fc-dependent inverse transverse gluon-propagator Ga 
obeys GA(g^) = GA(g^) +m\. 



3 Flow equation 

Our method amounts to a solution of flow-equations for the fc-dependence of the gluon and ghost 

two-point functions Ga(q'^) and Ggh{q^)^ starting at some high momentum scale k = ko with 

perturbative "initial conditions" and extrapolating to k —>■ 0. The exact flow-equations for the 

two-point functions obtain by taking the second functional derivative of equation (|]) and therefore 

involve the exact momentum dependent three- and four-point functions between gluons, quarks 

and ghosts. Those are not known and we have to proceed to some approximation. (In principle, 

the general flow equation (|]) also yields exact flow equations for the three and four point functions. 

Since the latter involve again higher n-point functions a solution of such equations would postpone 

the necessity of approximations to the higher n-point functions.) One possible approach would be 

the extraction of the three and four-point functions from the truncation (^, (^, expressing them 

in terms of the functions^ Zp, Yp and H . There is, however, no reason to believe that such a 

procedure is particularly accurate - it is easy to construct additional invariants contributing to the 

three and four-point functions as, for example, {P^^F'^^^ . Instead of such a truncation we will 
^See for the Hmit where Yp is neglected. 
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base our approximation on the observation that only a relatively narrow range in the integration 
momentum contributes effectively to the integrals on the right hand sides of the flow equations 
for dtGAiq"^) and dtGgh{q^)- For a given propagator momentum there should exist some effective 
gauge coupling such that a classical approximation of the renormalized vertices gives the same 
result as the full momentum dependent three and four-point functions. (This statement is close in 
spirit to the central value theorem for ordinary integration.) Of course, this effective coupling may 
be different for different values of and between dtGA^q^) and dtGgh{q^). Our approximation 
consists now in taking independent of the momentum of the propagators and using the same 
value for dtGA^q^) and dtGgh{q^). In practice this permits the use of the classical momentum 
dependence for the three and four-point functions on the right hand side of the flow equations. 
We should point out, however, that this does by no means imply the assumption that the vertices 
are close to the classical ones. We only suppose that a certain momentum weighted average can be 
represented by a similar average over classical vertices with appropriate gk- It is a drawback of this 
approach that in practice gt is only defined implicitly]], without an explicit representation in terms 
of n-point vertices. This results in an uncertainty about what should be the scale dependence 
of gk and we will specify our assumptions below. In fact, we believe that the main uncertainty 
concerns this scale dependence. The approximation of an effective g^-independence of gt probably 
plays only a minor role, at least for momenta above the confinement scale. 

In summary, we insert on the right hand side of the flow equations a classical momentum de- 
pendence for the three point functions and a constant four point function in terms of renormalized 
gauge fields = Z^a' A and ghosts = Z^j^^, = Z^j^^. The strength of the renormalized ver- 
tices is then determined by a common renormalized gauge coupling gk appearing in the covariant 
derivative — igkT^A^j^^. Our approximation can be compared with the so-called Mandelstam- 
approximation in the SDE approach with the important difference that in contrast to the SDE 

only a narrow momentum range contributes in the loop integrals involved in the flow equations. 

''It is not difficult, however, to write down a formal expression in terms of momentum weighted integrals over 
exact three and four point functions. 
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An effective "momentum averaging of vertices" should therefore work much better in our case. 

With our approximations the flow-equations for the two-point functions simplify considerably 
and one finds for Landau gauge-fixing a = 



(2vr) 



1 ^ 



V 



-^Zl{GA{p)+rn\ + Rk{p)) \GA{q+p)+rn\ + Rk{q+p)) 



X 



13g2 - lA{qp) + I2p^ - 12 



{qpf 

g2 



22 



{qpf {qpf g2 ^2^2 _ (^^^y 



1 ^. 



+ -ZI^Za {GgK{p) + Ru{p)r' {G,K{q+p)+Ru{q+p)r' 



q2p2 p2 (^g_|_p~)2 

[qpY' 



p 



(13) 



l-GAq') 



3glZAZ',^J-^M{GA{p)+m\ + Rk{p)y\G,f,{q+p) + Rk^^^ 



-1 , 



(27r 



2 _ (gP)^ 

pZ 



(14) 



Here the symbol dt denotes a t-derivative which acts only on the infrared cutoff R^, i.e. dt = 
{dRk / dt){d / dRk) ■ The exponential decay of dtRk (cf. equation (|T2|) ) for large p^/A;^ implies a fast 
ultraviolet convergence of the momentum integrals in (|13]) and ([1^). 

It remains to specify the scale-dependence of the coupling gk- Even though defined only 
implicitly, gk is a possible choice of a renormalized gauge coupling (in an appropriate implicitly 
defined scheme). For small gk its running must therefore be governed by the universal two-loop 
beta-function and we use for sufficiently large k 

d 



-22 



at 

167r2 



204 



9l 
(167r 



2^2 



for k > kr 



(15) 



(No quark loops are included here since the fermionic cutoff kp remains fixed. The running is 
then given by QCD with only heavy quarks or an effective number of light flavors Nf = 0.) For 
larger gk we have in the present approximations no knowledge about the evolution of gk- Also our 
approximation of a classical momentum dependence for the vertices becomes less accurate. We 
find it suitable to display the effect of these uncertainties by using a set of different /5-functions 
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for large gt- We replace ([T5|) for k < kchy 

id = c{gl -gl) ; c = ' 9c = alike) for k<K (16) 

5'c 5'* 

This has the effect that the increase of gk is smoothly stopped for small k. For k gk approaches 
a fixed value g^. The scale where the running of gk starts deviating from the two loop running is 
fixed by the condition g"^ = 0.95gl. Equation ([T6|) corresponds then to a one parameter family of 
/3-functions. We will parameterize the uncertainty of our approximation by 

with typical values of of order one. It will become apparent that our lack of knowledge of the 
exact vertices plays an important role for the shape of Viq"^) only at small q^. For ^^^600 MeV 
the dependence of l^(g^) on turns out to be weak. 

After performing the p integrations, the flow equations ([T3| ) and ( p!4D become a coupled system 
of nonlinear partial differential equations for two functions Ga, Ggh, depending on two variables 
and t. They are supplemented by equations (|^), (|T6D for the running of gk- We have solved these 
equations numerically, starting at ko from perturbative initial conditions for G'yi(g^) and Ggh{q^). 

4 Initial Values 

The "initial values" for the propagators at the scale fco are determined by the solution of the flow 
equations for k > ko, starting from very large k with the "classical action" including appropriate 
counterterms Tct- Since for k > k^ perturbation theory is a valid expansion we can treat the flow 
of Ga and Ggh for k > ko perturbatively. The flow-equations in the one-loop approximation are 
obtained by setting on the right-hand sides of ([T3|) and ([T^) G'a(q'^) = Ggh{q^) = q^, Zi = 1, 
m\ = and gl = g^. In this case the momentum-integrations can be performed explicitly. We 



find 



dGAiq^) fk'^\27 44 15 7x^^ , . ^ , , s/^^ / x Eii(x/2)~ 
- + ^ + — (Eii(x)-Eii(x/2))-a;MEii(x) 



dt IGtt^ [4 XX 

, / 9 9x 11 23 9 \ - / o 55 6 N I , 

' 1 1 X X^ J \ X x^ 
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and 



^ - ' +- + -a;2(Eii(x)-Eii(a;/2))+exp(-x) -— + - + - + 



dt 167r2 12x2''' ^ / ^ V 2 2 x 

-) 

X J 



+ exp(-x/2) ( 3a; - 6 - — U + 0{t) (19) 



where x = /k'^ and 



Eii(x) = ^ 



Note that the g^-derivatives of the above expressions reproduce for ^ the perturbative 
anomalous dimensions as they should. 

The coupling does not depend on k and we can integrate equations ( |l8l) and ([19D . This 
yields aX k = ko the one-loop expression^ (x = /k^) 

- + 77^ i TT (In X + Eii(x/2) ) + — + — - — - — (Eii(x)-Eii(x/2)) + 



q'^ IGvr^ [6 8x 3^^ 6x^ 12 

5 1 \ 

+ 



x^ / / X Eii(x/2) \ / X 2 3 5 1 \ 

l'i_l_±_ii4._LV 



and 

^^^^ = .,. + -|^||(lnx + Ei,(x/2)) + A+ 1 _^(Ei,(a;)-Eii(x/2)) + 
q^ loTT^ [4 4x 2x^ 4 

The integration constants za and z^/j are related to the choice of a normalization scheme (see next 
section) . For the gluon-mass counterterm we find to one loop |1^, ||] 



rn\ = ---^k^ (22) 
^ 8 167r2 ^ ^ 



We have also used equations (|^) - (|22|) as a powerful check of our numerical integration algorithms 

by solving numerically equations (p^, (|13]), and (|^) for small values of gl- 

^The integration of equations (|l^) and (jl^) should be performed between fcg and some ultraviolet scale Aq ^ 

2 

feg. For Ao/fco oo the contribution from the upper boundary Aq reduces to a simple term cx Eii(i^r2-) — s- 
— (lux + In . This yields the logarithms in equations ( [20| ) and (|2|) and a fco-dependent piece in the constants, 
ZA = const. — In ^ and similar for Zgh- 
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5 Renormalization schemes and the relation to the MS- 
scheme 



In order to give numbers in physical units we have to specify the renormahzation-scheme. For our 
purpose this specifies how a given numerical value for the dimensionless coupling g^^ is related to a 
scale ko in GeV units. The numerical integration of the fiow equations then yields all dimensionful 
quantities (as etc.) in units of ko. As emphasized in the introduction, the right hand side of 
the fiow-equation (H) is both IR- and UV-finite for a suitable choice of the cutoff-function i?^ such 
as the one used in this work. Thus the formalism of exact fiow-equations constitutes effectively a 
regularization-scheme. In this context the renormalization-scheme is defined by fixing a starting 
value for the functional at a large scale k = k^. Physical quantities involving typical momenta 
<^ ^0 should only depend on a few relevant (or marginal) parameters in F^g (up to corrections 
proportional to powers of q^ /k'^). In this sense the solution of the exact fiow equations must become 
independent of the renormalization scheme as ^ 0. On the other hand, truncated versions of 
the fiow equations will typically exhibit some dependence on the renormalization scheme. 

In addition to the relation between and k^ a complete determination of the renormaliza- 
tion scheme requires in our approximation also a specification of the initial values za and Zgh in 
equations (|20| ) and (0). We choose a parametrization 



167r2 I 24 2 24 ^ ^ ^ 12 | IQn^ ^ ' 



and 



z,H = l-^|--^ + Eii(l/2) + le-i-2e-i/4-V^ (24) 

such that GA,ko{q^ = k^) = Ggh,ko{,q^ = k^) = k^ for 5a = Sgh = 0. In principle, the dependence 
of the final results on details of the relation between g^.^ and ko or on 6 can be used as a possible 
test of the truncation. For k ^ the solution of the exact fiow equations should not show any 
scheme dependence for the physical quantities. A strong dependence of the approximative results 
on "scheme parameters" like 6 would therefore indicate a shortcoming of the truncation. (This is 
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quite similar to the scheme dependence in a given order in a loop expansion.) On the other hand, 
the choice of an "optimal renormalization scheme" may permit to keep the truncation errors on a 
minimal level. 

The most important specification of the renormalization scheme concerns the value of at 
the initial scale ko, where we remind that gk is defined only implicitly. It seems to be a good 
procedure to compute some physical quantity related to a mass scale in our scheme and to choose 
QkQ such that this quantity coincides with observation or the calculation in some other known 
scheme as the MS-scheme. We only want to use perturbative physics as an input and would like 
a simple quantitative relation between our scheme and the MS-scheme. A convenient physical 
quantity is the value of the Fourier transform of the effective heavy quark potential at a given 
momentum scale V{q'^ = /i^). The latter has been computed in the MS-scheme to two-loop order 



TH] and one has 



+ 



4343 , 97r^ 

^ + 5477^- — + 66C(3)- 

/1229 52^,^,\ /lO,, 
(^ + yC(3))iV,+ (-iV, 



(167r2)2 



(25) 



The four quark interaction V{q ) may differ from V{q ) by terms of order g /q (see section 6) 

V{q') = V{q') + ^ = V{q') + ^ (26) 
with A = A + A, cf. equation (^. We therefore determine Qk^ by imposing 



for a fixed value of fi in the perturbative range. 

We still need the value of A{fi^). To lowest order it can be inferred by comparing a one loop 
calculation of Vlq"^) in the MS-scheme with eq. (pS]). In one loop order the difference ^ — \^ is 
ultraviolet finite in Landau gauge. One obtains 



167r2 



A(/i2) = -5^^M|1^ (28) 
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and the relations (|27|) and ( p8|) are now sufficient to express g'^ = as a function of 5 and g'^{^). 
The nontrivial part of the solutions of flow equations concerns then the shape of V^(q'^) or V{q^) 
for 7^ /i^. 

Let us next turn to the choice of an optimal renormalization scheme which minimizes the 
importance of neglected corrections. For the choice of g (or /i in equation (p7|)) we rely on the fact 
that V is two loop exact for = /i^. All truncation errors in order g^ therefore only affect the 
behaviour of V for q"^ ^ fi'^. Since we are mainly interested in low values of we will choose /x^ 
as low as possible with the constraint that the two loop approximation V^''"°^ should still remain 
valid. As a typical scale we will take fi = 3A(^ with A^ the two loop confinement scale in the 
MS-scheme {Nj = 0). For an optimal choice of 6 we will further assume (see section 6) that the 
difference A(g^) — A(/i^) remains small as compared to q'^V{q'^) for a certain momentum range 

around fi^. This will allow us to relate the momentum dependence of Vlq"^) and Vlq"^). By 
expanding equations @, ( ^0|) and (^TJ) for large in powers of g"^ it is easy to convince oneself 
that the momentum dependence of g^V^(g^) is influenced in two loop order by the value of 6. 
With our approximations this also extends to q'^V{q'^). We will adopt a choice of 6 for which the 
deviation between V{q'^) and V^'°°^(g^) remains small for a certain range g^ > /i^. (Of course, 
a truncation which is exact in two loop order would fix the two loop momentum dependence 
unambiguously and V could deviate from V^'°°^ only in higher loop order. Our optimal choice of 
6 partly anticipates these higher order effects by minimizing their role in a whole range g^ > /i^.) 
For the practical calculations with Nf = we will take 6a = Sgh = 6 = —3. In the following 
we will always work with this "optimal renormalization scheme" . This turns the disadvantage of 
the scheme ambiguity into an advantage by making use of independent perturbative calculations 
for g^ > /i^. Of course, this optimal scheme does not tell us what will happen for g^ < /i^ and 
it is in this region where the solution of the flow equation should extend our information beyond 
perturbation theory. 

The choice of the scheme also affects the size of the neglected corrections ~ A(g^) — A(/i^) in 
the relation between V and V. The fact that this quantity vanishes identically in one loop order 
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and is almost zero for a whole range > fi"^ strongly suggests that this quantity is indeed small 
compared to q'^V{q'^) in the optimal renormalization scheme, at least for a certain range < 
We thus expect V{q'^) + to be a good approximation to the heavy quark potential in this 
range. On the other hand, the information on Vlq"^) can be transfered to a precise determination 
of the four quark interaction V{q'^) only once the term A in equation is computed. For the 
time being we identify V with V (i.e. we assume A = — j^). 



\ 


1 








































2-loop 














6 = -3 

























v^[GeV] 



Figure 3: Comparison of q'^V(q'^) and 5'^V^'°°^(g^) for different values of 5 and fi (see equations 
(22), (23) and (26)). The dotted line gives the two loop MS-result, the solid lines correspond to 
6 = —3 and the dashed lines are obtained for 6 = 4. The scale /x is choosen to be 3A|^ for the 
lower (upper) solid (dashed) line and 4.5A^ for the upper (lower) solid (dashed) line. 



Results for q^V{q'^) with two values /ii = 3A^ and fi2 = 4.5Af^ are shown in figure 3. (For 
the choice of the momentum scale see section 6). Since the present truncation of the renormal- 
ization group equations reproduces the running of q'^V{q'^) only to one loop accuracy, the scheme 
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dependence can still be relatively important, as can be seen by comparing our results for the two 
values of fi but with the same choice of S = 4. For a fixed value of fi we should, following the 
discussion above, find the value of 6 such that the two loop result is approximately reproduced 
in a range > /x^. Varying /i within this range then does not change the results considerably 
anymore, as can be seen from the curves with 6 = —3. 



6 Effective heavy quark potential 

It is tempting to associate V{q'^) — j^r;^ with the Fourier transform of the heavy quark potential 
V{q^). (See for an identification of V(g^) with Vlq"^).) The potential can then be inferred 
from the gluon and ghost propagators for A; as 



V 



Vi,')-, ^ ^11 .„.. -T^ (29) 



As we have argued in the preceeding section, we expect this to be a valid approximation in 
Landau gauge and for the "optimal renormalization scheme" used here. On the other hand, this 
identification is certainly not valid for general gauges and renormalization schemes. It is important 
to keep in mind the structural differences between the four quark interaction V and the heavy 
quark potential V . First of all, we observe that is a gauge invariant quantity whereas V and V 
are not. For example, the momentum dependence of q'^V{q^) coincides with the one of q^V{q^) in 
one loop order only in Landau gauge. 

The static potential is commonly defined in terms of Wegner- Wilson-Loops |]T3. An expansion 



of the expectation value of the Wegner- Wilson loop (or a similar gauge invariant quantity involving 
four quarks) in terms of Green's functions for quarks and gluons not only contains a lowest order 
four-quark interaction (||) but also higher terms oc g < {ipipyA >, g"^ < {ipipyA'^ > etc.. After 
solving the field equations for A such terms would appear in the form of effective interactions 
involving six or more quarks. They contribute to — in higher orders in g. 

Alternatively, one may think of using the effective four quark interaction for a computation of 
"heavy quark scattering" in the form of some physical and therefore gauge invariant cross section. 
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However, any reasonable definition of heavy quark scattering involves (almost) on shell quarks and 
must include soft or collinear gluons. In contrast, Vlq"^) is evaluated for Euclidean quark momenta 
and nonvanishing infrared cutoff kp- The relative suppression of quark- vertex corrections (Y^), 
box diagrams (E) and quark wave function effects {Z^ — 1) for low is due to the infrared cutoff 
kp- It may not be the same for on-shell massive quarks with kp = 0. (Only for Euclidean 
momenta we are guaranteed that quark propagators always give suppression factors oc kp^ or mjj^ 
with rriH the heavy quark mass if the diagram is ultraviolet finite.) 

In summary, we conclude that on the conceptual level the difference between V and V 
is far from being trivial. Nevertheless, we know that for Landau gauge fixing the difference 
q'^ (V {q'^) — V {q'^)^ = A(g^) can only contain terms oc g^, g^lnq^/jj?' etc. in a perturbative 
expansion. The constant term oc g'^ is incorporated in our approximation (^). The optimal 
renormalization scheme minimizes the combined error from the neglection of the g^-dependence 
of A and the truncation in the computation of V . Only in this scheme the neglected corrections 
to V remain small for a whole range q^ > fi"^, and, by continuity, also for a certain range < yU^. 
It seems then reasonable to propose that the solution of the flow equation yields in the present 
setting a reasonable picture of V(g^) down to relatively low values of g^, say ~ 300 MeV. 

In order to test the approximations we will compare our results for V^(g^) with phenomenolog- 
ical quark potentials used for charmonium, as for example the Richardson potential |^ 



Vniq') = ^- J -V , Aij = 0.4GeV (30) 



g21n(l + f 



Of course, such a comparison would be more meaningful if the effects of light quark flavors were 
included. Nevertheless, results of a "quenched approximation" {Nf = 0) should not come out too 
far from reality provided the overall mass scale is properly chosen. 

In this section we will present the results of a numerical integration of the system of integro- 
differential equations (0) and (|14D with a gluon-mass counterterm given by (|l^) and a running 
coupling as discussed in section 3. Since this amounts to a calculation without light quarks 
{Nf = 0) we compare our results to perturbation theory in a pure Yang-Mills theory. All mass 
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and momentum scales are expressed in terms of the two loop confinement scale Ajj^ for Nf = 0. 
The value of in physical units will depend on the problem under consideration. We concentrate 
in this section on a computation of the heavy quark potential in the quenched approximation. In 
order to relate to the real world where Nf = 3 we use the condition that the MS gauge 
coupling for Nf = should coincide with the one for Nf = 3 at the edge of the range of validity 
of perturbation theory. (This minimizes trivial effects from the different running of the coupling 
with and without light quarks in the perturbative range.) In practice, we identify the Nf = 
gauge coupling with the "real" gauge coupling at a scale where Pg^/g"^ = —1.5, resulting in 
^Ms ~ '^^^ MeV (for a three-flavour confinement scale A^~^'* = 285 MeV). For a computation 
of the effective quark interaction discussed in the introduction the choice of Aj^ will be different. 
As discussed in the conclusions, the Nf = gauge coupling should then coincide with the "real" 
Nf = 3 coupling at the scale k^. 

We start the numerical integration of the flow equations at ko = 20Aj^ = 6.3 GeV and use /i = 
3Aj^ = 945 GeV and 6 = —3 for the definition of the renormalization scheme in equation (pT]). For 
the infrared value of the gauge coupling we take = 1.5. We will first discuss the results for the 
heavy quark potential V(g^) as defined in (pQ]). In figure 2 we have displayed q'^V{q'^) as a function 
of \/q^ for various values of the infrared cutoff k. For k = 500 MeV the effects of fluctuations are 
still moderate since the low momentum fluctuations are suppressed. The value of q'^V{q^) remains 
near the classical value ^g'^{fi) ~ 11. As k is lowered to 100 MeV, additional fluctuations are 
included. One observes a strong rise of q'^V{q'^) for small g^. Finally, for — > all fluctuations are 
included and the solid curve gives our estimate for the effective potential. It is this latter curve 
which has been plotted for different values of fi and 6 and compared with perturbation theory in 
figure 3. For nonzero k we note that Vfe(g^) gives a good approximation to Volq"^) for ~ 3/c. 
Otherwise speaking, the infrared cutoff substantially influences the potential only in the region 
y/q^ ^ 3k. This is a more quantitative version of our statement that fluctuations with momenta 
q^^k"^ have only little influence on the n-point functions at large q^. In figure 4 we compare the final 
value of V{q^) with the Richardson potential for three values = 1, 1.5 and 2. The agreement 
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for — 1 - 1.5 is quite remarkable. Overall, we find that our results smoothly interpolate 
between the two loop potential for large and phenomenologically acceptable potentials for 
300 MeV < < 800 MeV. 
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Figure 4: Heavy quark potential for different values of a^,. The solid line corresponds to «^ = 1.5 
whereas the dashed lines are for a* = 1 (lower dashed line) and a* = 2 (upper dashed line). Also 
shown are the two loop potential and the Richardson fit as dotted lines. 



For a more detailed discussion of the momentum dependence of V{q'^) it is useful to define the 

quantity 

. d 



Xv{q 



In [q^V{q^ 



(31) 



din (f' 

We have plotted this quantity in figure 5, again for different values of A;. In the limit of a constant 

Xv one has 

A 



Xv^Vv 



V{q') 
24 



(32) 




Figure 5: Momentum dependent anomalous dimension Xvi'f) ^^i different values of the cutoff k. 
The dotted line corresponds to k = 500 MeV, the dashed line is at /c = 100 MeV, and the solid 
line displays the final result for A; = 0. 



A potential V{q^) oc which is the Fourier transform of a linear rising confinement potential 
corresponds to Xv{.<f 0) = —2. Our results have a tendency to approach this value. For 
V^^400 MeV the dependence on the truncation (i.e. the value of in equation (|1^) becomes 
important, however. It is apparent that the truncation uncertainty affects a derivative quantity 
as xv much stronger than the potential itself. In figure 6 we display the function Xv{(f ) for /c = 
for different values of a*. 

We observe that xv is directly related to the /5-function for a suitably defined running coupling 

d 

9v = Pgl{9v) = Xvgv (33) 



din \fq 
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Figure 6: Xvil"^), for /c = and different values of a^: tlie solid line is for a^, = 1.5, tlie dasfied 
lines are for a* = 1 (upper dashed line) and a* = 2 (lower dashed line). The dotted lines give the 
two loop value and the result corresponding to the Richardson potential. 



The two-loop perturbative value is given by the universal two-loop /^-function and thus reads 



X,,(,') = -22«IM- 20401 (34) 



Extracting gyiq"^) from our numerical solution we have also plotted the function ( P^ in figure 
6. The agreement of our numerical solution for xv with this curve for high is remarkable. 
We emphasize that the two loop running of in ([15]) does not automatically imply a two loop 
running of the momentum dependent gv{<f) for large q^. At the present level of accuracy it is 
impossible to distinguish which part of the deviation of Xv{q^) from (0) is due to shortcomings of 
the truncation or to the difference between the four quark interaction V{q^) and the heavy quark 
potential. 
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For small momenta, the perturbative result for xv is no longer reliable. For the Richardson 
potential (equation (0)), the behaviour of xv reads 

This curve is also displayed in figure 6. Again, we observe the interpolation between the two loop 
and Richardson potential. 

The quantity gy may also be used in order to parametrize our results for V{q^) in a form which 
allows for simple rescaling of the overall mass unit. Let us introduce the dimensionless variable 
z = g^/A|jg. Our results for the four quark interaction can be summarized in the fit (a^, = 1.5) 



47r^ ^ An ] 3 IGn'^ 



4343 o 9n 
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+ 5An' - — + 66C(3) 
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(167r2)' 



X 



f _ 3.283 7.262 3.584 1 , , 

x<^ 0.977+ - + ^ ^} 36 

[ 1 + 1.154Z (1 + 48.24^)^ (1 + 0.269^)^ J 



where 



ASn' j 102 In (ln(2: + 3.507)) I 1 
^^^^ 33 \ 121 ln(z + 3.507) J ln(;2 + 3.507) ^ ^ 

This fit reproduces our results in the range 0.05 ^ z ^ 100 with a relative error < 0.5%. 

Since the agreement with phenomenological potentials is satisfactory in the range 300 MeV < 
i/g^ < 800 MeV the charmonium levels computed from our approximation to Vlq"^) will be close 
to the experimental values. One may ask if a string tension a can also be extracted. Instead of 
determining a from the behaviour of Vlq"^) for — > we will use here a more phenomenological 
definition by fitting the quantity cr(g^) = ■^q^V{q'^) for ^A-^ < y/cp < |Aj^ as 

^iq') = ^q'V{q') = a + cq' (38) 
167r 

In figure 7 we display cr(g2) for different values of a^. For a range = 1 — 2 we extract 

0.13 GeV^ <a < 0.46 GeV^ (39) 
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Figure 7: String tension: The plot shows o{(^) = -^q^V{q'^) for different values of ay^. The solid 
line is for a^, — 1.5, the dotted line is for — 1 and the dashed line gives the result for a^, — 2. 
We also indicate the straight hne extrapolations from the region 0.025 GeV^ < < 0.06 GeV^. 



and 1.7 > c > 0.03. In view of the uncertainties and missing effects of light quark fluctuations the 
agreement with the observations from quarkonia (cr f« 0.18 GeV^) seems reasonable. Obviously, 
our approximations are insufficient to extrapolate V{(f) to — > and to establish a "confinement 
behaviour" V ~ unambiguously. This is apparent from the qualitative change of q^V{(f) in 
dependence on a* (figure 7). We observe a bifurcation in the sense that for > of' the limit 
of q^V{q^) diverges for ^ whereas for ct* < of* it vanishes. Only for a^, — af^ ~ 2 the 
potential scales exactly oc One may speculate that a proper treatment of the flow of the 

effective coupling could lead to a combined fixed point for gk (corresponding to a^''), ^t(?^ = 0) 
and Yp[q^ — 0), whereas Zp{q^ — 0) vanishes for /c ^ 0. It would be interesting to investigate 
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Figure 8: Wave function renormalizations Za (solid line) and Zgh (dashed line) as functions of 
ln(/c/Aj^), for a* = 1.5. 

the scaling behaviour of the flow equation for the ansatz (^ with Zf[c^^ = Kpq^ and Yp, Kp 
independent of q^. The vertices should then be defined by this truncation and the effective coupling 
Qk could be extracted a posteriori. 

Finally, it is also worthwhile to study the momentum and scale dependence of the gluon and 
ghost propagators in Landau gauge. The overall size of Ga and Ggh is given by the wave function 
renormalizations Za and Zgh (equations (|ll|))- The fc-dependence of Z{k) is plotted in figure 8 
and we observe in particular a strong decrease of Za for A; — > 0. The momentum dependence of 
the propagators is shown in figure 9. Our results seem to indicate that GA{q^)/q^ vanishes for 



0. 
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Figure 9: GA{q^)/q^ (lower lines) and Ggh{q^)/q^ (upper lines), for A; = and different values of 
a*: the solid lines are for = 1.5, the dotted lines correspond to a* = 1 and the dashed lines 
give the results for = 2. 



For a more detailed investigation of the momentum dependence we introduce "momentum 
dependent anomalous dimensions" 

, V M ^_2^ In '''''''' -^'<°' (40) 

The quantity xv in equation (^) is related to xa and Xgh by 

Xv{q') = {xA{q') + 2xM)) - ^^^^ (41) 

where we observe that the correction oc A becomes rapidly insignificant for small q^. We have 
plotted Xa and Xgh in figure 10. For large these quantities should approach the one loop 
anomalous dimensions for the wave function renormalization of the gluon and ghost fields. 
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Figure 10: Momentum dependent anomalous dimensions Xa(?^) (lower lines) and Xgh{(f) (upper 
lines), for A; = and different values of a*: the solid lines are for a* = 1.5, the dotted hues are for 
a* = 1 and the dashed lines for a^, — 2. 
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We have checked that our numerical solution obeys these relations. 



7 Conclusions 

We have computed the effective four quark interaction which results from QCD after integrating 
out the gluon (and ghost) degrees of freedom. More precisely, we consider an effective action where 
also the quark fluctuations with momenta larger than a fermionic infrared cutoff k^p are integrated 
out, but not those with lower momenta. In Landau gauge the low momentum behaviour of the 
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four quark interaction is dominated by a dressed one gluon exchange whereas contributions from 
box diagrams and the quark wave function renormahzation are less singular. In fact, the latter 
can be treated perturbatively since infrared complications are avoided due to the cutoff in the 
quark propagators. By the use of gauge symmetry (Ward identities) we can express the dominant 
"one gluon exchange contribution" V{q'^) in terms of the gluon and ghost propagators GA{q'^) 
and Ggh{q^)- These propagators have been computed by an approximate solution of an exact 
renormahzation group equation in the context of the effective average action. 

For high enough momenta there is a simple relation between V{q^) and the effective heavy 
quark potential 1^(5'^) in the quenched aproximation [Nj = 0). The difference can be computed 
perturbatively and we have taken into account the one loop contribution. Using this connection for 
Vq^ > 3Aj^ we have optimized our renormahzation scheme by employing the existing information 
on the two loop heavy quark potential in the MS scheme. This also allows us to express all 
dimensionful quantities in units of the two loop confinement scale Aj^. We should point out that 
no information on the behaviour of the heavy quark potential for low momenta < 3Aj;jg is 
required for this optimization. Our final result for V is the fit in equation (|36|) , 



The full four quark interaction V^(g^) (cf. equation (|)) still needs to be related to V^(g^) by a 
perturbative calculation, the difference being oc g^/q"^ and subleading for small g^. 

Our result (|43|) is expressed in terms of the dimensionless variable z = q'^/Ay^. In contrast 
to the discusssion of the quenched heavy quark potential in section 6, the effective four quark 
interaction depends on the scale k^. In fact, the running of the gauge coupling corresponds 
effectively to a flow with three light quarks, Nf = 3, for k > k^ and Nf = for k < k^. The 
low momentum range is dominated by the flow for k < k^ and the main dependence on k^ arises 
through the scale Aj^ in the fit (|36|). In fact, Aj^ increases with k^ due to the faster running of 
the gauge coupling for Nf = as compared to Nf = 3. (To one loop order one finds oc k"^, 
7 = 2/11.) For k^ = 1.3(1.7) GeV one obtains numerically Aj^ = 375(395) MeV. All momentum 
scales in the figures have to be rescaled correspondingly by a factor Ai^/315 MeV. With the 
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appropriate choice of our result (|43|) can now be used as an input for the flow of the effective 
quark interactions as the fermionic cutoff is lowered For future computations of the flow in 
the effective quark model an important check of the approximations consists in monitoring if the 
dependence on cancels out of the final results. 

We have also made an attempt to use our information on the effective four quark interaction 
K(g^) for a computation of the effective heavy quark potential V{q'^) in the quenched approxima- 
tion. This involves an additional important assumption, namely that the potential is dominated 
by a dressed one gluon exchange (in Landau gauge) also for ^/q^ < SAj^. In other words, we 
suggest that the difference between V^(g^) and Vlq"^) remains subleading for small q^. Even though 
suggested strongly by continuity arguments, this remains an assumption as long as no detailed cal- 
culations are available. In the real world the three light quarks contribute to the potential. (There 
is no infrared cutoff fc^ for these quarks.) In order to adapt our quenched results to this situation 
one has to choose the relevant scale appropriatly so that at least the high momentum fluctu- 
ations of the quarks are effectively included (see section 6). Comparing our results to observation 
(coded in phenomenological potentials) gives then very satisfactory results: We find an effective 
potential which interpolates smoothly between the two-loop perturbative potential for high mo- 
menta and phenomenologically successful potentials for the range 300 MeV^^/q^^800 MeV which 
is relevant for quarkonia. This strongly supports our assumption that Vlq"^) is indeed dominated 
by a dressed one gluon exchange. If this picture can be substantiated and extends to the straight- 
forward inclusion of light quark fluctuations it will permit to relate the perturbative scale Aj^ to 
quarkonium physics and therefore ultimately compute the fine structure constant at the Z-mass 
as{Mz) in the MS scheme. Although the error estimate is difficult and probably not competitive 
with a corresponding program in lattice QCD such an analytical link between perturbative and 
nonperturbative QCD would be of great value. 
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